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The kinetic origin of resonance phenomena in capacitively coupled radio frequency plasmas is discov¬ 
ered based on particle-based numerical simulations. The analysis of the spatio-temporal distributions 
of plasma parameters such as the densities of hot and cold electrons, as well as the conduction and 
displacement currents reveals the mechanism of the formation of multiple electron beams during 
sheath expansion. The interplay between highly energetic beam electrons and low energetic bulk 
electrons is identified as the physical origin of the excitation of harmonics in the current. 


Capacitively coupled radio frequency (CORE) discharges 
are indispensable tools for semiconductor manufactur¬ 
ing and other innovative applications [U |2]. At the same 
time, they are challenging physical systems due to their 
complex and nonlinear dynamics. At low neutral gas 
pressures of a few Pa or less, CORE discharges are op¬ 
erated in a strongly non-local regime. In the so-called 
“(a-mode”, electron heating is dominated by stochas¬ 
tic sheath expansion heating [3] and electric field re¬ 
versal during sheath collapse [IHZl- Stochastic heating 
was modelled extensively in the past in the frame of 
a hard wall model, as well as pressure heating EHni. 
During the phase of sheath expansion, energetic electron 
beams are generated and propagate into the plasma bulk, 
where they sustain the discharge via ionization and lead 
to a Bi-Maxwellian electron energy distribution function 
(EEDE) [T8ti24] . 

At low pressures, resonance effects such as the plasma 
series resonance (PSR) [27II3Q] and the plasma parallel 
resonance (PPR) j3TH33] can be self-excited and strongly 
enhance the electron heating mug. In the presence of a 
sinusoidal driving voltage waveform, the excitation of the 
PSR results in a non-sinusoidal RE current, due to the 
appearance of harmonics of the driving frequency [30] . 
Although these have been observed in experiments ED, 
they are usually neglected in most models of electron 
heating in CORE plasmas. Existing theories which in¬ 
clude resonance effects are zero-dimensional global mod¬ 
els based on equivalent electrical circuits [29] as well as 
spatially resolved models based on the cold plasma ap¬ 
proximation m- As these models do not include any 
kinetic effects, such resonances should be investigated on 
a microscopic kinetic level. A kinetic interpretation is re¬ 
quired to clarify some of the most important open ques¬ 
tions about electron heating dynamics in CORE plas¬ 
mas: What is the kinetic origin of the generation of high 
frequency (HE) oscillations of the RE current and the 
generation of multiple electron beams during one phase 


of sheath expansion such as observed in previous works 
[37] ? In what way is current continuity (V • jtot = 0) 
ensured at all times within the RE period in the pres¬ 
ence of electron beams, where the total current density 
itot = jd + jc is decomposed into the displacement and 
conduction current density? 

Our aim is to provide access to a kinetic interpreta¬ 
tion of resonance phenomena in CORE discharges based 
on Particle-in-Cell simulations complemented by Monte- 
Carlo treatment of collision processes (PIC/MCC). We 
consider a single frequency (/ = 55 MHz) low pres¬ 
sure {p = 1.3 Pa) argon plasma driven capacitively by 
a sinusoidal voltage source. The simulation code is 1- 
dimensional in space and 3-dimensional in velocity space. 
The electrodes are assumed to be infinite, planar, paral¬ 
lel, and separated by a gap of L = 1.5 cm. One of the elec¬ 
trodes (at X = 0) is driven by a sinusoidal voltage wave¬ 
form, 0(t) = 00 sin (2 co’rf^) 5 where cjrf = 27r-55 MHz and 
the other electrode {x = L) is grounded. A symmetrical 
setup is adopted to demonstrate that resonance phenom¬ 
ena are of general relevance even in symmetric discharges, 
where they are usually neglected. Eirst, we analyze the 
discharge behavior at 0o = 150 V driving voltage ampli¬ 
tude and subsequently, we discuss the discharge behavior 
at 00 = 300 V. Secondary electron emission and particle 
reflection at the electrodes are neglected in order to sim¬ 
plify the interpretation of the electron heating dynamics. 
The cross sections for electron-atom and ion-atom colli¬ 
sions are taken from [35j [36] . Eor more details see [37] . 

Eigure shows the applied sinusoidal driving voltage 
and the calculated current waveform that is found to be 
nearly 90 degrees phase shifted and non-sinusoidal as a 
consequence of the self-excitation of the PSR [T^[2^[3Q] . 
The vertical dashed lines indicate characteristic reference 
times which are used in the following discussion. We 
note that due to their inertia (cjpi <C cjrf, where 
is the ion plasma frequency) the motion of ions is not 
modulated with the excitation frequency. Ions follow the 
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FIG. 1. Driving voltage waveform 0o = 150 V (left scale) and 
calculated current waveform (right scale) for one RF period. 
The vertical dashed lines indicate reference times (ti, t 2 , ts, G) 
used in the forthcoming analysis. 


time-averaged electric field. 
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FIG. 2. Top panels: Density of hot electrons {s > 11 eV) 
that move upwards (a) and downwards (b). Bottom panels: 
Density of cold electrons below 4 eV that move downwards (c) 
and upwards (d). The time axis covers one RF period. Den¬ 
sities are provided in m“^. The vertical dashed lines indicate 
the reference times specified in Figure 1. The powered elec¬ 
trode (150 V) is at the bottom (at x = 0), while the grounded 
electrode is at the top (at x = L). 


In order to unveil the kinetic origin of the non-linear 


current waveform, it is necessary to separate the dynam¬ 
ics of cold and hot electrons. The energy threshold for 
hot electrons is taken to be 5h = 11 eV, which is in the 
vicinity of the first excitation level of Ar atoms, while 
cold electrons are defined to have an energy below 5c = 
4 eV. The energy is calculated using all three velocity 
components {sh,c = Figureshows 

the density of electrons partitioned according to their 
energy (as explained above) and according to their direc¬ 
tion of velocity. Electrons with Vx > 0 (i.e., those moving 
from the powered electrode towards the grounded elec¬ 
trode) are defined to move “upwards”, while electrons 
with Vx < 0 are defined to move “downwards”, in ac¬ 
cordance with the representation of the spatio-temporal 
distributions displayed in Figure The first row of Fig¬ 
ure shows spatio-temporal plots of the number den¬ 
sities of hot electrons with energies above 11 eV only, 
which move upwards [Fig. (a)] and downwards [Fig. 
|^(b)]. Similarly to previous works [HI [37], the genera¬ 
tion of multiple beams of energetic electrons is observed 
adjacent to each electrode during the phase of sheath ex¬ 
pansion [here: 2 pronounced and 1 weak beam. Figs. 

(a) and (b)]. The bottom row presents the density of cold 
electrons with energies below 4 eV. These results are in 
strong contrast to the prevailing picture of the genera¬ 
tion of a single electron beam during one sheath expan¬ 
sion phase. The generation of multiple electron beams 
strongly affects the spatio-temporal ionization and exci¬ 
tation dynamics, which show similar structures m (i-e., 
this phenomenon is essential for the generation of such 
plasmas). Figure [^includes the same vertical lines as Fig¬ 
ure indicating the same characteristic times within one 
RF period. Figurej^shows the momentary electron veloc¬ 
ity distribution function (EVDF) at these times, spatially 
averaged over the plasma bulk [6 mm < x < 9 mm, Eigs. 
|^(a) - (d)]. In each plot, electrons with energies above 11 
eV are represented by red bars, electrons with energies 
below 4 eV are marked in green and electrons with ener¬ 
gies between 4 eV and 11 eV are marked in blue. Here, 
the energy is calculated using only the component (ax¬ 
ial direction), since the Vy and Vz components do not cou¬ 
ple due to the fact that electrons are accelerated only in 
axial direction by the expanding sheaths and almost no 
collisions cause an angular scattering of electron beams. 
The dashed lines in Eig. correspond to a Maxwellian 
distribution fitted to the data shown in Eig. (a). 

In the following, a detailed kinetic interpretation of the 
electron heating dynamics during the phase of sheath 
expansion at the bottom electrode is presented. (Due 
to the symmetry, the phenomena occur at the top elec¬ 
trode half an RE period later.) During sheath collapse 
(ti ^ 4.4 ns) at the bottom electrode the total current is 
zero and the EVDE is approximately Maxwellian. Eol- 
lowing the start of the sheath expansion (^2 ~ 8.1 ns), the 
first beam of energetic electrons is generated and propa¬ 
gates towards the plasma bulk [Eig. [^. These beam elec- 
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FIG. 3. The momentary electron velocity distribution func¬ 
tion (EVDF) at the considered reference times, spatially aver¬ 
aged over the plasma bulk (6 mm < x < 9 mm). The dashed 
lines indicate a Maxwellian distribution, which is htted for 
t — ti. Electrons above 11 eV are represented by red bars, 
electrons below 4 eV by green bars and electrons between 4 
and 11 eV by blue bars (00 = 150 V). 


trons move away from the expanding sheath edge leading 
to a strongly anisotropic EVDF [Fig. (^ 2 )]- Clearly, 
the number of electrons at high positive velocities is in¬ 
creased. As the beam electrons move away from the ex¬ 
panding sheath edge, they leave a positive space charge 
behind. This positive space charge causes an electric field 
that accelerates electrons back towards the sheath edge. 
Consequently, shortly after (ts « 10 ns) cold bulk elec¬ 
trons move back towards the sheath edge [Fig i (^3)]. 
At this time, energetic beam electrons move upwards, 
while cold bulk electrons move downwards [see Fig. 

(c)] (i.e., two groups of electrons move into opposite di¬ 
rections simultaneously). The bulk electrons cannot re¬ 
spond instantaneously to the perturbation caused by the 
energetic beam electrons due to their inertia. They can 
only respond on the timescale of the local electron plasma 
frequency (cjpe ~ Scjrf in the bulk), which is approxi¬ 
mately 3.5 ns (r ~ 27r/ci;pe). This leads to a modulation 
of the electron density in the bulk at the local electron 
plasma frequency. (Recall that the inertia of electrons 
is the kinetic reason why an “inductance” is required 
in global equivalent circuit models to excite resonance 
effects.) Later (^4 « 11.6 ns), the drifting cold bulk 
electrons approach the energy barrier of the expanding 


sheath edge and upon impact a second beam of ener¬ 
getic electrons is generated [Fig. [^(^ 4 )]. This process is 
repeated until the sheath expansion stops. Under these 
conditions, this mechanism leads to the generation of two 
pronounced energetic electron beams during one phase 
of sheath expansion at t ~ 8.1 ns, 11.6 ns. Furthermore, 
a third weak beam formation at 15 ns is present. For 
the latter the modulated bulk electrons reach the end 
of the sheath expansion and, therefore, experience just 
a slight kick by the sheath potential. This kinetic pic¬ 
ture links the complex dynamics of the electrons to the 
non-linearity of the RF current. 
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FIG. 4. Spatio-temporal distribution of the conduction (a) 
and displacement (b) current density for one RF-cycle in 
A/m^ (00 = 150 V). 



FIG. 5. The conduction and displacement current density 
(left) for one RF-cycle and the corresponding Fourier spectra 
(right) in the center, at 7.5 mm (0o = 150 V). 

Figure shows spatio-temporal plots of the conduc¬ 
tion and displacement current density. The presence of 
energetic beam electrons causes a local enhancement of 
the conduction current density at distinct times within 
the RF period [Fig. |^(a)]. In order to ensure current 
continuity (V • jtot = 0 ) the plasma must react to this 
local perturbation instantaneously. Consequently, a dis¬ 
placement current is generated, where the beam prop¬ 
agates [Fig. (b)]. This local displacement current is 
180° out of phase with respect to the conduction current 
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and compensates the perturbation caused by the pres¬ 
ence of the beam electrons locally [Fig. |^. This can be 
understood as a local plasma parallel resonance (PPR) 
[3TI433] . Physically, this displacement current is gener¬ 
ated by the beam electrons themselves, since they move 
away from the expanding sheath edge. They leave be¬ 
hind a positive space charge, which the bulk electrons 
cannot instantaneously compensate. This space charge 
in turn causes an electric field, that first increases and 
then decreases, when the space charge is compensated 
by the cold bulk electrons. In this way current continu¬ 
ity is ensured in the presence of ballistic beam electrons 
based on this kinetic mechanism. In the voltage driven 
case, the displacement current compensates the conduc¬ 
tion current only partially at higher odd harmonics of 
the driving frequency [Fig. [^. Thus, certain harmon¬ 
ics in the total current, and consequently, the PSR, are 
self-excited. (Even harmonics are not excited due to the 
symmetry of the discharge.) 

Increasing the driving voltage amplitude (0o = 300 V) 
leads to a higher central plasma density and electron 
plasma frequency (cjpe ~ Ocjrf in the bulk). The cold 
bulk electrons can respond faster to the perturbation 
caused by the energetic beam electrons [Fig. [^. They 
move back towards the expanding sheath edge on shorter 
timescales (r ~ 2 ns) and the above process can be re¬ 
peated more often during one phase of sheath expansion. 
Consequently, the number of electron beams generated 
during one phase of sheath expansion increases [Fig. 
(a)] at higher driving voltage amplitudes. In the latter 
case, the kinetic interplay between low and high energetic 
electrons leads to the excitation of higher harmonics (9th 
harmonic) in the total current [Fig. as well as in the 
current densities in the center of the discharge, which can 
be seen from the waveform and spectra given in Fig. 

In summary, the mechanisms discussed here correspond 
to the first kinetic interpretation of resonance phenom¬ 
ena in CCRF plasmas. It is demonstrated that these 
effects cannot be understood in detail neither based on 
fluid models nor on global equivalent circuit models, since 
(i) local phenomena play a crucial role and (ii) a kinetic 
analysis is required to describe two groups of electrons 
propagating in opposite directions simultaneously. In 
principle, these kinetic mechanisms are present in any 
CCRF discharge (asymmetric or symmetric) and they are 
key processes that drive the electron heating dynamics. 
However, their experimental identification, particularly 
in the given symmetric case, is challenging due to the 
high temporal resolution required (on a sub-nanosecond 
timescale). 
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FIG. 6. Density of hot electrons {e > 11 eV) that move up¬ 
wards (a) and density of cold electrons below 4 eV that move 
downwards (b)^ 
shown in Fig. 


(b) in m ^ (0o = 300 V). Similar color scale as 
5 . Dfrom 0 to a) 4.05 10^^ and b) 1.21 10^^ m” ~ 



150 


100 


50 


-50 


-100 


-150 


FIG. 7. Driving voltage waveform 00 = 300 V (left scale) and 
calculated current waveform (right scale) for one RF period. 



FIG. 8. The conduction and displacement current density 
(left) for one RF-cycle and the corresponding Fourier spectra 
(right) in the center, at 7.5 mm (0o = 300 V). 
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